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ABSTRACT 

We report the discovery of a galaxy cluster at z = 1.41. ISCS J143809+341419 was found in 
the Spitzer/IRAC Shallow Survey of the Bootes field in the NOAO Deep Wide-Field Survey carried 
out by IRAC. The cluster candidate was initially identified as a high density region of objects with 
photometric redshifts in the range 1.3 < z < 1.5. Optical spectroscopy of a limited number of 
objects in the region shows that 5 galaxies within a ^120 arcsec diameter region lie at z = 1.41 ± 
0.01. Most of these member galaxies have broad-band colors consistent with the expected spectral 
energy distribution of a passively-evolving elliptical galaxy formed at high redshift. The redshift 
of ISCS J143809+341419 is the highest currently known for a spectroscopically-confirmed cluster of 
galaxies. 

Subject headings: galaxies: clusters — galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

High-rcdshift galaxy clusters provide important 
tools in the study of galaxy formation and evolu- 
tion. The galaxy populations of rich cluster cores 
at z > 1 tend to be dominated by massive ellip- 
ticals (e.g. iStanford et~aTl Il997t iRosati et al l | 1999| 
Stanford et all |2002i Ivan Dokkum fc Stanford! 120031 
BlakeskTe e^ al. ||2Q03j: |Rosafo et^ aL||2004: Postma n et alJ 
20051 iMullis et all 120051: iMei et alJ 120051) . which are 
useful probes of galaxy evolution because their stellar 
populations appear to be relatively simple. As the 
highest overdensity regions, clusters should contain 
the oldest (and most massive) galaxies, so that by 
probing z > 1 we can determine useful constraints 
on the galaxy formation process. The connection 
between such stellar population s and the ellipti- 
cals in present-epoch clusters ijBower et alJ 119921) 
has been the subj ect of extensive analysis in re- 
cent years (e.g. | ^a nford 1 Eisenharcjt^ fc, Dickinscinl 



T993 Ivan Dokkum et al.l 119981 



De Prop ris et alJl!999t iKelson et al 



_ubi n et all [1998; 
2000]). The nature 
of elliptical galaxy formation in clusters at z > 1 de- 
pends strongly on the importanc e and mode of merging 
in assembling the ste llar mass ijvan Dokkum fc Franxl 
120011 lConselice"ll2005D . Beyond z ~ 1 in cosmologically- 
flat CDM models, the amount of merging occurring 
within the prior ~1 Gyr is large and should seriously 
inflate the locus of early-type galaxy colors du e to 
interaction-induced starbursts (Kauffmann 1996]), al- 
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though evidence ex ists for "dry" mer gers that could 
alleviate this issue (jvan Dokkuml l2005j) . On the other 
hand, models in which ellipticals formed by a monolithic 
collapse at high-z predict a tight color-magnitude 
relation out to at least z ~ 2 for reaso nable cosmologies 
ijEggen. Lvnden-Bell. fc Sandagelll9"62Th The identifica- 
tion of clusters at z > 1 and the characterization of their 
galaxy populations provides a powerful means of testing 
elliptical galaxy formation theories. 

One proven method of finding z > 1 cluster candidates 
is deep infrared sky surveys. Because they are massive 
and their rest-frame light peaks at 1 — 2/im, early-type 
galaxies stand out from the background in the observed- 
frame near- and mid-IR light out to z ~ 2. Therefore, to 
the extent that galaxy clusters are composed of galaxies, 
near- and mid-IR imaging surveys offer a viable alter- 
native to X-ray searches for very high redshift clusters. 
Both methods need to be exploited to sample the entire 
range of galaxy clusters in the key 1 < z < 2 epoch be- 
cause it is unclear when and how massive galaxies and 
the associated intracluster medium (ICM) are formed. 
One advantage of using infrared imaging surveys is that, 
in addition to identification of massive clusters with an 
ICM that would be detected by an X-ray survey, they can 
be extended to probe to lower mass scales and identify 
the building blocks that yield massive clusters at lower 
redshift. We describe here the discovery of a z > 1 cluster 
which resulted from an imaging s urvey conducted by the 
InfraRed Array Camera (IRAC; iFazio et aT] 120041) on - 
board the Spitzer Space Telescope l|Werner et alJ l20041 . 
Except where noted, the assumed cosmological parame- 
ters are H = 70 km s" 1 Mpc" 1 , O m = 0.3 and A = 0.7. 

2. OBSERVATIONS 

2.1. Sample Definition 

To identify galaxy cluster candidates in the Bootes 
field we made use of the IRAC Shallow Survey 
ijEisenhardt et alJl2004 - the FLAMEX Survey, which 
provides JK S imaging down to K s — 19.2 (R. Elston 
et al. in preparation), and the NOAO Deep Wide-Field 
Survey (NDWFS), w hich provides deep B W RI imaging 
ijJanuzzi fc Devlll999|) . A catalog was defined using the 
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Fig. 1.— I,K,4.5 fim false-color image of ISCS J143809+341419 
covering 3' on a side. Objects with photometric redshifts at 1.25 < 
z < 1.55 are marked by the yellow boxes. Photometric redshifts 
were not calculated for some of the fainter objects in the 4.5/tm 
catalog so some objects visible in the image which are not marked 
with boxes in fact may be in the cluster. The 5 spectroscopically 
confirmed member galaxies are marked by the red circles. North 
is up and East is to the left in the image. 



4.5/im data from the IRAC survey. Photometric red- 
shifts were estimated from matched photometry in the 
B W RIJK and the IRAC 3.6 and 4.5 /jm bands. Then 
a wavelet search was carried out using the photometric 
rcdshift probability distributions to detect structures on 
physical scales of 200 < r < 800 kpc. Bootstrap simula- 
tions were then used to quantify the significance of each 
detection. The details of the sample selection and iden- 
tification of cluster candidates are given in P. Eisenhardt 
et al. (in preparation), and of the photometric rcdshift 
process in M. Brodwin et al. (in preparation). We report 
here on the candidate cluster ISCS J143809+341419, 
shown in Figure ^ The prefix ISCS stands for the IRAC 
Shallow survey Cluster Search. The cluster overdensity 
relative to that of similar size areas in the entire Bootes 
field in the IRAC Shallow Survey is shown in Figure [21 

2.2. Keck Spectroscopy 

In order to make sure that the candidate cluster is not 
a chance alignment along the line of sight, we obtained 
optical spectroscopy of objects in ISCS J143809+341419. 
Objects were selected based on their photometric red- 
shifts. The /-band magnitudes of the primary galaxy 
targets were 22.0 < I < 24.3. We prepared a slitmask 
including slitlets for 14 objects with photometric red- 
shifts 1.25 < z p < 1.55 within 3 arcmin of the nominal 
cluster center. The range in z p was chosen to cover the 
nominal photometric rcdshift of the cluster Ha. Two 
other objects in the 4.5/im catalog were included so as 
to fill out the mask. The slits had widths of 1.3 arc- 
sec and minimum lengths of 10 arcsec. The slitmask 
was used with the Low Resolution Imaging Spectrograph 



Fig. 2. — Histogram vs rcdshift showing the photometric red- 
shifts of galaxies in the IRAC Shallow Survey. The area covered 
by the solid red histogram is 3' X 3' around the z = 1.41 cluster. 
The solid black histogram represents the average of 100 3' X 3' 
patches randomly chosen from the entire ~9 degree 2 IRAC survey 
in the Bootes field. The dashed histograms represent the 1 and 2 
a uncertainties on the average value histogram. 



ijOke et al.lll995|) on the 10 m Keck II telescope on UT 
2005 June 04 to obtain deep spectroscopy. On the red 
side, we used the 400 1 mm" 1 grating, which is blazed 
at 8500 A, to cover a nominal wavelength range of 6000 
to 9800 A, depending on the position of a slit in the 
mask. The dispersion of ~1.8 A pixel -1 resulted in a 
spectral resolution of R ~ 900. On the blue side we used 
the 400 line grism which is blazed at 3400 A and pro- 
vides coverage from the atmospheric cutoff up to ~5800 
A where the dichroic splits the light between the two 
sides of LRIS. We obtained 7 x 1800 s exposures with 
this setup in photometric conditions with 0.9 arcsec see- 
ing. Objects were shifted along the long axis of the slits 
between exposures to enable better sky subtraction and 
fringe correction. The observations were carried out with 
the slitlets aligned close to the parallactic angle. 

The slitmask data were separated into individual spec- 
tra and then reduced using standard longslit techniques. 
A relative flux calibration was obtained from longslit ob- 
servations of the standard stars HZ 44 and Wolf 1346. 
One-dimensional spectra were extracted from the sum 
of all the reduced data for each of the 16 slitlets for both 
the red and blue sides. For the targets in the cluster, 
only the red side data are useful and the blue side data 
will not be discussed further. 

3. RESULTS 

3.1. Optical Spectroscopy 

Spectra were obtained for 9 of the 14 photometrically- 
selected galaxies with quality sufficient for determining 
the redshifts; the other 5 were deemed too faint for reli- 
able identification of spectral features. Only one of these 
9 objects has a redshift outside of the photometric red- 
shift range. Four objects have redshifts within the range 
1.4147 < z s < 1.4172 so are considered to be mem- 
bers of a cluster. A fifth object has z s — 1.4028 which 
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Fig. 3. — Histogram of all redshifts determined in the Keck/LRIS 
slitmask observation. The shaded objects are those with photomet- 



is somewhat lower than the redshifts of the other four 
but still likely to be a cluster member. A histogram 
of all redshifts determined from the mask (including 5 
serendipitous sources) is shown in Figure El The spec- 
troscopy validates our photometric redshifts in the range 
1.3 < z < 1.5 (M. Brodwin et al. in preparation). 

The spectroscopic members are listed in Table 1. The 
table gives the Vega magnitude [4.5] in the IRAC 4.5 /im 
band, the photometric and spectroscopic redshifts, and 
the template number of the best-fit as determined by the 
photom etric redshift est imator. These template numbers 
ar e fromlBrodwinl (l2004f) based on th e spectral templates 
of iColeman. Wu. fc Wccdman (1980). The approximate 
correspondence to galaxy type spectra is for elliptical, 
7 for an Sbc, and 13 for an Scd. In line with the tem- 
plate fitting, the optical spectra of most of the members 
show absorption lines (Ca II H+K, A3830, Mgl A2852 
and Mgll A2800) and spectral breaks (D4000, B3260, 
B2900) similar to those of present-epoch ellipticals. The 
spectra of the member galaxies are shown in Figure 
Only one of the 5 members,IRAC J143811+341256, has 
a prominent emission line; the photometry of this same 
object was fit best by a relatively late-type template. To 
better determine the redshifts of the 4 early-type clus- 
ter members, we used the Fourier quotient technique 
as implemented in the FXCOR task of IRAF. The red 
members were cross-correlated with an early-type galaxy 
template (Kin nev et al.lll99fi[) . The resulting redshifts, 
along with their formal uncertainties, are listed in Ta- 
ble 1. We adopt a nominal center for the cluster of 
a = 14 ft 38 m 09?4, 5 = +34°14'19" (J2000), the cen- 
troid of the wavelet detection, and a nominal redshift 
z = 1.413. Given the small number of members it is pre- 
mature to estimate the velocity dispersion of the galaxies. 

3.2. Optical and Infrared Photometry 

A color-magnitude diagram for all objects in a 3' area 
at the cluster are shown in Figure [S] The objects with 
photometric redshifts in the range 1.25 < z < 1.55 are 
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Fig. 4. — Optical spectra of the 5 member galaxies after be- 
ing smoothed by an 5 pixel boxcar. The following spectral features 
(though not present in every object) are marked in increasing wave- 
length by the dotted lines: Mg IIA2800, B2900, B3260, [OII]A3727, 
Mg IA3830, Ca II K and H, and D4000. The rest frame wavelength 
at z = 1.4166 is shown along the top. 
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Fig. 5. — Color-magnitude diagram of the cluster in / — [4.5] 
vs [4.5]. The magnitudes are in the Vega system and [4.5] is the 
magnitude in the 4.5 ^m band. The spectroscopically identified 
galaxies at z = 1.41 are marked with blue squares, and the photo- 
metrically selected galaxies at 1.25 < z < 1.55 by red diamonds. 
The green horizontal line represents a model prediction at L* for 
a passively evolving stellar population formed in a single 100 Myr 
burst beginning at Zf = 3.0, as described in the text. The errors 
in the photometry are shown only for the spectroscopic member 
galaxies. The photometry represents 5" corrected aperture magni- 
tudes. The field size is 3'x 3'. 



shown by the red diamonds, and those with spectroscopic 
redshifts 1.40 < z < 1.42 are shown by the blue squares. 
The scatter seen in the colors of the member galaxies is 
dominated by the photometric uncertainties in the ND- 
WFS photometry of I > 23 galaxies. 

Also plotted in Figure is an estimate of the ex- 
pected color for L* early-type galaxies. Using a simple 
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passive evolution model c alculated from the GISSEL of 
iBruzual fc Charlotl l)2003l) and our assumed cosmology, 
we calculated the expected colors for a single 100 Myr 
burst stellar population formed at Zf — 3.0. As can be 
seen in FigureEl the colors of the brighter member galax- 
ies are consistent with those exected for early-type galax- 
ies whose stars formed at Zf — 3.0. The three brightest 
objects on the color-magnitude relation have not yet been 
observed spectroscopically. 

4. DISCUSSION 

The observations presented here provide 
strong evidence supporting the identification of 
ISCS J143809+341419 as a galaxy cluster at z = 1.41. 
Photometric redshifts, demonstrated to be accurate by 
the optical spectroscopy in the relevant redshift range, 
indicate a significant overdensity of galaxies within a 
region of radius r — 300 kpc at z ~ 1.4. Our optical 
spectroscopy shows that 5 of these galaxies are at the 
same redshift. Further examination of the properties of 
the cluster, such as its mass and the evolutionary state 
of its galaxies, is deferred until spectroscopy can confirm 
more members and deeper multiband imaging can be 
obtained. 

It is noteworthy that ISCS J143809+341419 was de- 
tected using imaging taken by IRAC with exposures last- 
ing only 90 s per pointing over the ~9 degree 2 area in 
the Bootes field of the NDWFS. As implied by the iden- 
tification of another z ~ 1.4 cluster whic h relied on a 
short 20 ks exposure in the XMM archive ijMullis et all 
2005), current and planned surveys are better off cov- 
ering larger areas rather than going deep if the goal 
is to find z > 1 galaxy clusters. To date three other 
z > 1 IRAC-selected clusters have been spectroscopi- 
cally confirmed in the Bootes field (P. Eisenhardt et al. 
in preparation, M. Brodwin et al. in preparation, and R. 



Elston et al. in preparation). With the combination of 
the search described here in the Bootes field, serendip- 
itous searches for clusters in the XM M archival data 
ijRomer et al.ll2001tlMullis et al.ll2005D. and the on-goin g 
wide-area surveys such a s SWIRE ( Lonsd ale et all2 0031 , 
the RCS (|Gladders fc Yedl2005ft . and the CFHT-Legacy 
Survey, the time finally is ripe for the identification of 
large samples of z > 1 clusters. The construction of such 
samples will pave the way towards a better understand- 
ing of the origin of early- type galaxies. 
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TABLE 1 

Summary of Spectroscopic Cluster Members 



ID 


R.A. a 


Dec. a 


[4.5] 


phot-z 


template 


spec-z 


IRAC J143806+341433 


14:38:06.97 


34:14:33.8 


16.68 


1.39 


1 


1.4153 


IRAC J143807+341441 


14:38:07.80 


34:14:41.6 


16.18 


1.30 


3 


1.4166 


IRAC J143813+341452 


14:38:13.36 


34:14:52.9 


16.93 


1.48 





1.4028 


IRAC J143812+341318 


14:38:12.10 


34:13:18.2 


16.90 


1.35 


3 


1.4147 


IRAC J143811+341256 


14:38:11.25 


34:12:56.2 


17.20 


1.39 


11 


1.4172 



Coordinates arc J2000. 



